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The microscopic origin of superconductivity has 
been established in numerous classes of materi- 
als. In elemental metals it results from the ex- 
change of phonons whereas in copper oxides and 
iron pnictides it is intimately connected to mag- 
netism. On the other hand, the cause of super- 
conductivity in Bai a;Ka;Bi03 ( =32 K)[lJ and 
a large third class of isotropic high temperature 
superconductors is still mysterious and subject of 
debate. [2j Proposed mechanisms in bismuthates 
involve either the tendency of the Bi atoms avoid- 
ing the Bi^+ state, while preferring the Bi'^+ and 
Bi^+ charge state thus generating an attractive 
local "negative U" pairing center [3, 4j, or strong 
coupling to an O-breathing phonon mode[5l [6]. 
However, first principles calculations have so far 
ruled out both the negative U mechanism[7l [8] 
and the electron-phonon mechanism [9j. Here we 
demonstrate that dynamical correlations among 
electrons enhance their coupling to phonons in 
a large number of superconductors, raising sig- 
nificantly the previously underestimated theoret- 
ical superconducting critical temperature Tc of 
up to a few Kelvin to their experimental values 
of tens of Kelvin. The mechanism we propose 
is quite general, explains the magnitude and the 
doping dependence of superconductivity in many 
materials such as the celebrated Bai 2;Kj;Bi03 
compounds and the electron-doped /3-HfNCl com- 
pounds (Tc = 25.5 K)10 and can be used to design 
other high temperature superconductors which 
have not yet been synthesized. As an applica- 
tion, we suggest Bai^La3;Pb03 with x ^0.7 will 
superconduct above 10 K. 

Earlier approaches to the electron phonon coupling 
(EPC) were based on the density functional theory 
(DFT) and linear response theory jTT] in the local den- 
sity approximation (LDA) or generalized gradient ap- 
proximation (GGA). This approach correctly accounts 
for the rather high temperature superconductivity in ele- 
mental lithium (Tc -20 K)[T2], yttrium (T^ -20 K)[T3], 
calcium (Tc ~25 K)[T3] and binary compound MgB2 
(Tc ~40 K)fT5], so the results in ref [9 of rather small 
EPC strength A ~0.3 and very low Tc (up to 5 K) appar- 
ently rule out the phonon mechanism in Bai-^rKajBiOa. 

The widely used LDA/GGA exchange-correlation 



(EXC) functional in the DFT methodology, however, can 
suffer from large systematic errors in certain materials. 
These errors have been traced to two sources: failure 
to capture dynamic correlations because of self interac- 
tion errors or equivalently in the chemical sense failure to 
capture the delocalization of electrons due to Coulomb 
interaction (so-called delocalization errors in ref. jl6]). 
and failure to capture static correlation effects resulting 
from degeneracies. [TF Treating static correlation effects 
requires multi-determinant methods such as DFT com- 
bined with dynamic mean field theory, whereas the dy- 
namic correlations can be treated by using more sophis- 
ticated functionals such as hybrid functionals (HSE)[17] 
with nonlocal exchange potential. It has been recently 
pointed out[TSl [TS] that many difhculties in accounting 
for the physical properties of the BaBiOa related ma- 
terials within standard DFT method, are due to dy- 
namic correlations, and can be remedied by using hybrid 
functionals [ITj or GW calculations which largely remove 
the LDA/GGA self interaction (dynamic correlation) ef- 
fects. These improvements have resulted in excellent 
agreement with experiment for many computed physical 
properties, including crystal structure, lattice vibrations 
and optical conductivity. [El [19] 

Here we build on the insights of these works [T51 
[TO] to address the problem of the superconducting 
mechanism in these materials. Using DFT and lin- 
ear response theory]!!] with both the conventional 
LDA/GGA functional and the non-local hybrid func- 
tional (HSE06)[T7],(see supplementary material for cal- 
culational details) we evaluate the EPC, and find that 
the conventional functionals underestimates the EPC and 
predict the wrong trend with doping. On the other 
hand we find that once dynamic electronic correlations 
are taken into account, the EPC is sufhciently strong to 
account for the high temperature superconductivity in 
these materials. 

The middle atom M in BaBiOs (BBO) and BaPbOa 
(BPO) is octahedrally coordinated by 6 oxygen atoms as 
in cuprate superconductors. Ideally without any distor- 
tion (eg. when optimal doped such as Bao.6Ko.4Bi03), 
they crystallize in the simple cubic structure as shown 
in Fig. Jl] However, the octahedral formed by oxygen 
atoms deviates from the ideal octahedral in the parent 
compounds. [20l |2T] There are two types of distortions: 
one is the oxygen breathing distortion along nearest- 
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FIG. 1: The crystal structure of BaBiOa and BaPbOa. 

(a)The oxygen breathing mode/distortion. (b)The oxygen 
tilting mode/distortion. Both oxygen breathing and tilting 
distortions exist in BaBiOs compound whereas only oxygen 
tilting distortion exists in BaPbOa compound. 



neighbor M-0 bonds as shown in Fig. [T]^a). The other 
one is the oxygen tilting distortion perpendicular to the 
nearest-neighbor A/-0 bonds as shown in Fig. [l]^b). Ex- 
perimentally, both distortions exist in BBO whereas only 
0-tilting distortion exists in BPO. 
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FIG. 2: Band structure of BaBiOa and BaPbOs. The 

band structure of (a) BaBiOa and (b) BaPbOs calculated us- 
ing both the GGA and HSE06 hybrid functionals with the 
ideal undistorted simple cubic structure, the experimental 
structure, and for BaPbOa an O-breathing distorted struc- 
ture with O displacement of 0.085 A. Both GGA and HSE06 
produce similar band structure at the ideal structure around 
the Fermi level. For BaPbOa, the difference between GGA 
and HSE06 is less significant thus only GGA result is shown. 



The deformation potential D = Ae^/AQ is an inex- 
pensive way to measure the EPC strength, where Sk is 
the band eigenvalue at k point and Q is the vibration 
amplitude. Large deformation potential is an indica- 
tion of possible strong EPC. One example is provided 
by MgB2. In first principle calculations, the deformation 
potential, however, depends on the EXC functional, sig- 
naling the crucial role of EXC functional in evaluating 
the EPC strength in real materials. 

To demonstrate the EXC functional and Af-atom 
dependence of the deformation potential, we show in 
Figj2ja) the band structure of BaBiOa and BaPbOa with 
experimental and ideal structures using both GGA and 
HSE06 functionals in the folded face-centered cubic Bril- 
louin zone. The deformation potential of BaBiOs from 
HSE06 functional of the conduction band (mostly of Bi 
6s and O 2p states) across the Fermi level around L 
point, where the O-breathing instability occurs, is about 
13.3 eV/A, which is almost twice of the GGA value 7.6 
eV/A, and is even larger than the big deformation po- 
tential 13 eV/A of the B E2g breathing vibration mode 
found in the 40 K superconducting MgB2 compound. [25] 
In BaPbOa compound, however, the same deformation 
potential given by GGA and HSE06 bears smaller dif- 
ference, with values of 10.1 and 11.2 eV/A, respectively. 
The GGA deformation potential depends more strongly 
than HSE06 on the atom bonding to oxygen, specifically, 
the number of electrons of its outer shell. Using the vir- 
tual crystal approximation, we find that changing the 
charge of the M atom, the GGA deformation potential 
increases roughly linearly from A/=Bi to Pb, whereas 
changing the charge of Ba atom barely changes the de- 
formation potential. 

The mode EPC strength At, can be estimated through 
the deformation potential, following ref.j^ 

A. = ^(^.)^ E kv4P-i.4 (1) 

^ " j = l,2,3 

for HSE calculations of the O-breathing mode in BaBiOs, 
where the sum on j runs over the three moving O atoms 
of the O-breathing mode whose vibration frequency is 
ujv w65 meV, and iV(ei?)=0.36 eV/states is the density 
of states per spin at the Fermi level contributed by the 
corresponding band whose average deformation potential 
on the Fermi surface is Dj. 

Alternatively, we can estimate the EPC strength 
by comparing the deformation potential of HSE with 
GGA/LDA functional coupled with the A calculated from 
linear response calculations using GGA/LDA functional. 
Such an estimation is based on the fact that, the mode 
dependent EPC strength is proportional to squared de- 
formation potential, i.e., X-u oc D^, and the total A is 
usually contributed mostly by only a few phonon modes. 
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We therefore has a simple relation of 



A 



D 



HSE 



HSE 



'D 



,2 ^ 



LDA/GGA 



* -^LDA/GGA- 



(2) 



The superconductivity in doped BaBiOs is mainly con- 
tributed by the 0-breathing mode, analogous to the case 
of MgB2 where the superconductivity is mainly driven by 
the B-breathing vibrations. [22], The estimated Xhse ~ 
1.01 based on the above deformation potentials and 
Alda = 0.33 estimated from linear-response calculation, 
which is in agreement with previous LDA linear response 
calculations. [9] Using the modified McMillan equation 



T, 
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1.04(1 + A) , 
A-/i*(l + 0.62A)' 



(3) 



with experimental average frequency wjc 



450 K[ 



and A ~ 1.01, is estimated to be 31.8 K with /Lt*=0.1 
which is obtained by our first principle calculations (see 
supplementary material). As a result, the strong EPC 
strength A ~ 1.0 is enough to explain the rather high 
superconducting temperature Tc ~ 32 K in K-doped 
BaBiOa in the framework of the conventional phonon 
mediated superconductivity. 
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FIG. 3: Electron-phonon coupling strength A and su- 
perconducting Tc of Bai-^La^^PbO^. The doping effect 
is based on virtual crystal approximation. A (LDA/GGA) is 
calculated from linear response and AfHSE) is inferred from 
A(LDA/GGA) according to equation Q. 

Equation ([2| can be further generalized to two isostruc- 
tural compounds. As shown in Table |lj the relation 

Af^^/GGA - * ^f^A/GGA holds for the op- 

timal doped isostructural BaBiOa and BaPbOa. 

According to equation ([2|, the optimal Xhse for 
electron-doped BaPbOs is estimated to be 0.72, result- 
ing in a superconducting Tc=16.7 K, assuming u;;og=450 
K and /i*=0.1. Therefore, superconductivity (T^ up 
to 16.7 K) is expected in electron-doped BaPbOs com- 
pounds. This is already confirmed experimentally by 



doping Bi/Sb into the Pb site, resulting in Tc ~ 12 K for 
BaPbo.75Bio.25O3 [24] and 3.5 K for BaPbo.75Sbo.25O3 
compound[25j. Here we show in Figjsjthat superconduc- 
tivity may also be induced by doping into Ba site via 
studying the doping dependence of A in Bai_2:Laa;Pb03 
based on virtual crystal approximation. The strong dop- 
ing dependent A suggests it is crucial to bring the Fermi 
level close to the bands with large deformation potential 
in order to optimize superconducting T^. 

The mechanism of dynamical correlation (particularly 
in BaBiOs nonlocal exchange) enhanced phonon medi- 
ated superconductivity seems very relevant to many ma- 
terials where the electron bonding has both partially co- 
valent (or metallic) and partially ionic character. There 
are precisely the materials where high temperature su- 
perconductors have been found. Our observation of en- 
hanced EPC gives new insights into numerous mysteries. 

In heavy fermion and copper oxide based materials, 
two dimensionality enhances superconductivity [26] but 
in BaBi03 related materials, layering by synthesis of 
Ba„+iBi„03„+i {n~l, 2, 3, • • • ) degrades the critical 
temperature instead. We have evaluated the deforma- 
tion potential associated with 0-breathing distortion and 
found it to be almost zero for the Ba3Bi2 07 compound 
(n=2 member), resulting in very weak EPC and non- 
superconductivity, as seen in experiments |27j . For the 
Ba2Bi04 compounds (n—1 member), the corresponding 
deformation potential is about 70% of that in BaBiOs. 
Optimal doped Ba2Bi04 compounds could be supercon- 
ducting around 8 K, however, it hasn't been synthesized 
successfully in this form. 

The mechanism described in this letter, has a broad 
range of applicability. Relatively high temperature su- 
perconductivity was discovered in electron doped MNX 
(Af=Ti, Zr, and Hf; N=nitrogen; X=C1, Br, and I). T^ is 
16.5, 16 and 25.5 K for electron-doped a-TiNCl, /3-ZrNCl 
and ^-HfNCl, respectively. [351 [2S] ■ 

Also here, linear response calculation based on LDA 
functional produced A about 0.52 and ujiog about 36.4 
meV.[5n] which gives Tc only 6 K assuming ^*=0.1 and 
is insufficient to account for its 16 K superconductivity, 
raising the question of what is the mechanism of super- 
conductivity. A proper theory should explain why /3- 
HfNCl has higher Tc than /3-ZrNCl and why Tc in these 
electron-doped MNX compounds depends weakly on the 
doping level and the X atom. In other known high tem- 
perature superconductors, including the BaBi03 based 
ones, the critical temperature depends crucially on the 
level of doping. Our calculations indicate that in the 
electron-doped /3-ZrNCl compound, the deformation po- 
tential of the conduction bands (see supplementary ma- 
terial) associated with Zr-N vibration along c direction 
is substantially enhanced from 4.0 eV/A in GGA to 4.7 
eV/A in HSE functional. The enhanced deformation po- 
tential results in A ~ 0.72 and Tc=15.7 K (assuming 
a;/og=36.4 meV and fi*=0.1), in close agreement with ex- 
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TABLE I: Various quantities for selected compounds. Deformation potential D (eV/A) from GGA and HSE functionals 
for the most important vibration mode, electron- phonon coupling constant A from linear response calculations (LDA/GGA) 
and inferred value for HSE according to equation ([2|, average phonon frequency uiiog (K) from linear response calculations 
(LDA/GGA) and experiments, the calculated optimal superconducting temperature Tc (K) by LDA/GGA and HSE functionals 
(assuming /i*=0.10) and experimental Tc for selected compounds with proper electron or hole doping. 



Compounds 


mode 


D (GGA) 


D (HSE) 


A(LDA/GGA) 


A (HSE) 




Tc(LDA/GGA) 


Te (HSE) 


Tc (exp.) 


BaBiOs 


0-Bi-O 


7.6 


13.3 


0.33 


1.01 


450 


0.5 


31.8 


32.0 [I] 


BaPbOs 


O-Pb-O 


10.1 


11.2 


0.58 


0.72 


450 


9.3 


16.7 


12.0 [24] 


BasBiaOr 


0-Bi-O 


~0 


~0 


~0 


~0 


~450 


~0 


~0 


<2 [27] 


ZrNCl 


Zr-N 


4.0 


4.7 


0.52 


0.72 


422 


6.0 


15.7 


16 [H] 


HfNCl 


Hf-N 


4.4 


5.3 




0.91 


422 




25.2 


25.5 [To] 


TiNCl 


Ti-N 


3.8 


4.6 




0.70 


422 




14.7 


16.5 [21] 



perimental value of 16.0 K. More importantly, the same 
deformation potential is substantially larger in /3-HfNCl 
than in /3-ZrNCl, which is increased from 4.4 eV/A in 
GGA to 5.3 eV/1 in HSE. Since /3-H£NCl is isostructural 
to /3-ZrNCl, the 5.3 eV/A deformation potential gives A 
- 0.91 and Tc=25.2 K for electron-doped /3-HfNCl, in 
good agreement vi^ith experiment. It explains naturally 
the higher in ;3-HfNCl than /3-ZrNCl. The 16.5 K su- 
perconductivity in electron-doped a-TiNCl can also be 
well accounted for in the same way as shown in Table |l] 
Therefore linear response calculations with hybrid func- 
tional can fully account for the observed superconduct- 
ing Tc in these electron-doped MNX compounds. We 
expect that the unexpected high Tc superconductivity in 
many other superconductors including those reviewed in 
ref. (31 can be explained in a similar way by incorporat- 
ing appropriate EXC functional. 

As a final remark, we stress that first principle calcu- 
lations based on LDA/GGA functional can significantly 
underestimate the EPC strength due to the failure to 
capture dynamic correlations. As a result, critical reex- 
amination of the role of EPC in many superconductors 
including cuprate superconductors and recently discov- 
ered iron-based superconductors is called for by applying 
appropriate EXC functional such as HSE. The super- 
conductivity in many exotic superconductors therefore 
could be reconciled with phonon-mediated mechanism. 
Our study demonstrate that deformation potential can 
serve as a fast and effective criteria for searching potential 
phonon-mediated high Tc superconductors and doping is 
crucial to optimize Tc in this class of systems. 
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Method 

The electronic structures are calculated using both 
the VASP[gi] code with GGA(PBE)[£2] and HSE06[g3] 
exchange-correlation functionals and WIEN2k[S4 code 
with GGA(PBE)[^ exchange-correlation functional. 
Linear response calculations are carried out using the 
LMTART[35l |6] code with the ideal cubic structures 
(i.e., without distortions). The crystal structures are 
taken from ref.S[7] for BaBiOg and ref. S[H] for BaPbOg. 
For VASP calculations, the PAW-PBE pseudopotential 
is used, which is also coupled with the HSE calculations. 
The energy cutoff for the wavefunction is 400 eV. The re- 
sults of GGA calculations from VASP are double checked 
with results from WIEN2k with GGA(PBE) exchange- 
correlation functional. In linear response calculations, a 
32^ k mesh is used to converge the electron charge den- 
sity. The electron phonon coupling constant is obtained 
with a 8^ q mesh with 35 independent q points in the 
irreducible Brillouin zone. We use virtual crystal ap- 
proximation to simulate doping effect. 



The resulting n* for the ideal BaBiOa, BaPbOa and 
LaPbOs compounds are 0.104, 0.117 and 0.119 respec- 
tively. 

The difference between our work and Ref. 321 is that 
we use finite temperature formalism and approximate S 
function with the imaginary part of LDA Matsubara's 
Green's function 



i5(ekA - Ep) w Im- 



1 



(7) 



where ujq is the smallest positive Matsubara's frequency. 
The expression Q becomes exact in the limit of zero 
temperature. Technically we perform the following steps. 
First we calculate LDA one-electron spectrum and cor- 
responding wave functions. It provides us with LDA 
Green's function. Then we evaluate LDA polarizabil- 
ity Plda = GldaGlda- The third step is to calculate 
screened interaction W = V+VPldaW, where V is bare 
Coulomb interaction. After the third step we are able to 
get the Coulomb pseudopotential parameter. 



Evaluating Coulomb pseudopotential 

To estimate the Coulomb pseudopotential ji* , we im- 
plemented the methodology of Ref.[S9j in our codes [3T0l 
[TT| which combine LDA and GW methods in a linearized 
augmented plane- wave (LAPW) basis set. fj* is given by 
the formula 



l + /iln5 



(4) 



where e^^ is the Fermi energy and lod the Debye energy, 
and dimensionless parameter /i is calculated as the fol- 
lowing average over Fermi surface 



Ekk' Eaa' T^kA;k'A"5(ekA - ^^F)^(ek'A' ~ Ep) 



Ekk' Eaa' ^(f^kA - Ep)S{ei,'x' - Ep) 



(5) 



where N{Ep) is the density of states at Fermi level per 
spin, EkA is the Kohn-Sham eigenvalue of the Ath band 
at a wave vector k, and matrix elements of the screened 
Coulomb interaction in a basis of LDA eigenstates are 
calculated as the following 

W^kA;k'A' = (k'A' t, -k'A' ; |14^|kA t, -kA i). (6) 



Density of States 

The HSE06 functional also describe the doped com- 
pounds better. We show in Fig |Sl| the calculated den- 
sity of states (DOS) of an optimal doped compound 
Bao.62Ko.38Bi03 using both the PBE functional and 
HSE06 functional, based on virtual crystal and rigid band 
approximation, and we compare them with the experi- 
mental ultraviolet photoemission spectra [9T2] (UPS) of 
Bao.62Ko.38Bi03. We note that the rigid band approx- 
imation is as good as the virtual crystal approximation 
since we use the experimental structure. There is an over- 
all good agreement between the DOS from IISE06 calcu- 
lation and experimental UPS, in particular they share 
roughly the same peak centered at about -3.1 eV. In con- 
trast, the DOS from PBE calculation has a peak centered 
at about -2.2 eV. 



Optical conductivity 

The optical conductivity offer a baseline to assess the 
electronic structure and importance of electronic corre- 
lations. We plot in Fig. S2 the optical conductivity for 



the parent BaBiOs and BaPbOa compounds (with distor- 
tions) and the optical conductivity and effective number 
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FIG. SI: Density of states of Bao.62Ko.38Bi03. The DOS 

of Bao.62Ko.38Bi03 calculated using both the PBE functional 
and HSE06 hybrid functional. The experimental data is the 
ultraviolet photoemission spectra taken from Ref. 3121 



of charge carrier for the 0.38 K-doped Bao.62Ko.38Bi03 
compound in the cubic structure. 

For the parent BaBiOs compound, the DFT-PBE peak 
position in the optics shows a clear deviation (0.5 eV) 
from the experimental value as shown in Fig. |S2[ a). The 
HSE06 functional over-corrects the DFT-PBE peak po- 
sition, resulting in an overestimation of 0.3 eV. The best 
match is obtained by confining the non-local Fock ex- 
change within certain distance through a control param- 
eter /i=0.5, for which HSE functional is denoted as mHSE 
in Fig. |S2| Note that yit = corresponds to no confine- 
ment of the Fock exchange, i.e., the PBEO functional, 
fj, = oo corresponds to no Fock exchange is included, 
which restores to the PBE functional, while fi — 0.2 lies 
between the two limits and corresponds to the HSE06 
functional. Similarly, the HSE06 functional improves 
over DFT-PBE in the optics of BaPbOs compound as 



shown in Fig. S2 b) 



For the Bao.62Ko.38Bi03 compound, by comparing 
the experimental and DFT-PBE calculated optical con- 
ductivity and the energy-dependent effective number of 
charge carrier shown in Fig. |S2[ c) and (d), we estimate 
an overall mass enhancement of m* /niband ~ 2 in the 
low energy region, in consistent with the specific heat 
measurement of m* /mband up to 1.8.[S16j As discussed 
in the main article, the strong electron-phonon coupling 
(A ~1.0) in the optimal doped BaBi03 gives rise to a 
mass enhancement of m* /rriband = 1 -I- A ^ 2 thus is 
able to account for the above mass enhancement. There- 
fore, the effect of static electronic correlation in BaBi03 
is mild. 



o 
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FIG. S2: Optical conductivity and effective charge car- 
rier number, (a) and (b): the optical conductivity for the 
parent compound of BaBi03 and BaPbOs, respectively, (c) 
the optical conductivity and (d) the resulting effective carrier 
number Ng// for Bao.62Ko.38Bi03. We compare the calcu- 
lated results obtained by using both the PBE functional and 
HSE hybrid functional with available experimental data which 
are taken from Ref. 313H15I 



Lattice distortions and phonon softening 

We compute the total energy as a function of the two 
main distortions, i.e., 0-breathing and 0-tilting distor- 
tions. The results are shown in figure [S3] in the sto- 
ichiometric compounds. The inclusion of nonlocal ex- 
change has a dramatical effect on the lattice dynamics in 
BaBiOs but less so in BaPb03. Compared to standard 
PBE calculations, the 0-breathing instability is dramat- 
ically enhanced in the mHSE and HSE06 calculations 
in BaBi03 as suggested by the energy barrier [9T7] (0.25 
eV/f.u. in HSE06 vs almost eV/f.u. in PBE) shown 



in Fig S3 (a) and (b), whereas there is no 0-breathing 
instability in BaPb03 in both PBE and HSE06 calcula- 
tions as shown in Fig S3 (d) . For the tilting distortion 
on the other hand, the energy barriers (^ O.leV/ f.u., 
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FIG. S3: Total energy under distortions. The change of 
the total energy as a function of O-breathing distortion and 
O-tilting distortion from its ideal positions obtained by PBE 
functional and HSE hybrid functional (HSE06 and mHSE) in 
DFT calculations for BaBiOa compound (a, b, c) and BaPbOa 
compound (d). The distortions are normalized to its experi- 
mental distortion, i.e., /S.x / /S.Xexp- Since experimentally there 
is no O-breathing distortion in BaPbOs compounds, its O- 
breathing distortion is normalized to the experimental distor- 
tion in BaBiOa compound. 



Fig S3 (c) and (d)) are substantially smaller than the 
energy barrier 0.25 eV/f.u. of O-breathing distortion 
in HSE calculations. Note the PBE calculations have 
slightly larger energy barrier than HSE calculations for 
the tilting distortion. These results suggest that for the 
BaBiOs compound, the HSE functional strongly soft- 
ens the O-breathing phonon frequency compared to PBE 
functional. This is reflected in the average phonon fre- 



quency UJlog^ which in the LDA calculation uiog is 
about 550 K, while in experiments [3T9l [20] which re- 
semble the HSE functional calculation, is substantially 
smaller, only 450 K. 



Band structure of /3-HfNCl 

The DFT-PBE band structure and density of states 
of /3-HfNCl is shown in Fig|S4[ which is calculated with 
its experimental crystal structure without/with Hf 
atom shifted along the c-axis to change the Hf-N bonding 
length. In the calculations, the experimental stacking of 
the layers along the c-axis is ignored for simplicity. As the 
arrows indicated in the left panel of Fig |S4[ there is sub- 
stantial deformation potential for the first two conduc- 
tion bands, which are very flat along the F-Z symmetry 
line, resulting in a high density of states as shown in the 
right panel of FigjS4j Such a high density of states makes 
electron-doping hardly change the chemical potential in 
this region. The flat bands with substantial deforma- 




Z 1 2 3 4 5 
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FIG. S4: Band structure and density of states of /3- 

HfNCl. The band structure and density of states of /3-HfNCl 
calculated using the PBE functional with the experimental 
crystal structure [S2.1j and with the experimental crystal struc- 
ture except the Hf atom is shifted along the c-axis such that 
the Hf-N bond is shortened by 0.09 A. For simplicity, the 
different layer stacking along c-axis is ignored. 



tion potential is therefore likely the reason that the Tc 
is weakly doping-dependent in electron-doped /3-HfNCl 
compounds and other doped AfNX compounds. 
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